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Abstract 



Hadronic observables in Z+jet events can be subject to large NLO corrections 
at TeV scales, with X-factors that even reach values of order 50 in some cases. 
We develop a method, LoopSim, by which approximate NNLO predictions can be 
obtained for such observables, supplementing NLO Z+jet and NLO Z+2-jet results 
with a unitarity-based approximation for missing higher loop terms. We first test 
the method against known NNLO results for Drell-Yan lepton pt spectra. We then 
show our approximate NNLO results for the Z+jet observables. Finally we examine 
whether the LoopSim method can provide useful information even in cases without 
giant K-factors, with results for observables in dijet events that can be compared 
to early LHC data. 
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1 Introduction 

At CERN's Large Hadron Collider (LHC), it is widely anticipated that signals of new 
physics, for example supersymmetry, may manifest themselves as large excesses of data 
compared to expected QCD and electroweak backgrounds at high momentum scales PQEl 
El in 13 El E] ■ The estimation of these backgrounds will be one of the elements in ascertain- 
ing the presence of any new physics from such signals. Consequently, considerable effort 
is being invested across the particle physics community in the development of methods 
to understand and predict backgrounds (some of the issues involved are described nicely 
in ref. [5]). 

Given the QCD methods that are available today, some of the best prospects for 
obtaining systematic, accurate predictions of backgrounds involve next-to-leading order 
(NLO) QCD calculations. By carrying out a systematic expansion in the strong coupling 
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Figure 1: The LO and NLO distributions obtained with MCFM 5.7 [12] for three ob- 
servables in Z+jet production: the Z transverse momentum (left), the pt of the hardest 
jet (middle), and the scalar sum of the transverse momenta of all the jets, -ffTjets (right). 
The bands correspond to the uncertainty from a simultaneous variation of hr = fip i>J 
a factor of two either side of a default /x = Vpfji + m|. The jet algorithm is anti-fct [20] 
with R = 0.7 and only events whose hardest jet passes a cut pt > 200 GeV are accepted. 
The cross sections include the branching ratio Z — )■ e~^e~ . 



and obtaining the first two terms (leading order (LO) and NLO) for a given process, 
one often obtains predictions that are accurate to 10 — 20%. The importance of NLO 
predictions in the LHC programme has motivated a large calculational effort destined to 
extend the range of processes known at NLO (for reviews, see refs. [9l [TO]). 

While the majority of NLO calculations show some degree of convergence relative 
to the LO results, several groups have commented in recent years on the appearance of 
K factors, ratios of NLO to LO results, that grow dramatically towards high transverse 
momenta [TH [T^l [121 [HI [151 US] (similar behaviour is visible also in [T71[TS]). The problem 
generally occurs for hadronic observables (jet transverse momenta, etc.) in processes that 
involve heavy vector bosons or heavy quarks, at scales far above the boson or quark mass. 

Fig. [U illustrates this for the pp — )■ Z+jet process at LHC (14TeV) energies. It shows 
the distributions of three observables that are non-zero for configurations involving a 
Z-boson and one or more partons: the transverse- momentum of the Z-boson {pt,z), the 
transverse-momentum of the highest-pt jet {pt,ji) and the effective mass (scalar sum of 
the transverse momenta) of all jets (i^Tjcts)- At LO, all three distributions are identical. 
At NLO, the pt^z observable is rather typical of a QCD observable: its distribution has 
a NLO i^-factor of about 1.5, fairly independently of pt^z, and its scale dependence is 
reduced with respect to LO. The ptji distribution is more unusual: at high pt it has 
a ii'-factor that grows noticeably with Ptji, reaching values of about 4 — 6, which is 
anomalously large for a QCD correction. The -f^rjets observable is even more striking, 
with i^'-factors approaching 100. 

Given that fig. [1] involves momentum scales where as ~ 0.1, one is driven to ask how 
it is that such "giant" i^'-factors can arise. As touched on in pS], and discussed in more 
detail in [TH[T^ for the ptji case, the answer lies in the appearance of diagrams with new 
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Figure 2: A) a LO contribution to Z+jet production; B) and C) two contributions that 
are NLO corrections to Z+jet observables but whose topology is that of a dijet event with 
additional radiation of a soft or coUinear Z-boson either from a final-state quark (B) or 
an initial-state one (C). 

kinematic topologies at NLO. This is illustrated in fig. |2) at LO the only event topology 
(A) is that of a Z-boson recoiling against a quark or gluon jet. One type of NLO diagram 
involves gluon radiation from this basic topology, giving modest corrections to all our 
observables. However, there are also NLO diagrams (B,C) whose topology is that of a 
dijet event, in which a soft or coUinear Z-boson is radiated from outgoing or incoming 
legs. These diagrams do not contribute significantly to the pt^z distribution, because the 
Z-boson carries only a moderate fraction of the total pt- However when examining Pt,ji, 
it is irrelevant whether the Z boson is soft or not. Contributions B and C then lead to 
a result that is of order a^agw In^ptji/mz, where the double logarithm comes from the 
integration over soft and coUinear divergences for Z emission. The ratio of the NLO to 
LO results is therefore O (^agln'^ Pt,j±/mz)^ rather than just 0{as), hence the i^'-factor 
that grows large with increasing pto For the H^jets observable the enhancement is even 
bigger because the dijet topology leads to -f^Tjets ~ "^Ptji instead of -f^Tjets = Ptji at LO. 

While it is reassuring that we can understand the physical origins of the large K- 
factors in fig. [H we are still left with doubts as to the accuracy of the NLO Z+jet 
predictions for ptji and i^Tjets, since they are dominated by the LO result for the Z+2- 
parton topologies. One way forward would be to calculate the full NNLO corrections for 
the Z+jet process. However, while work is progressing on NNLO calculations of 2 — )■ 2 
processes with QCD final states (see e.g. |22j and references therein), results are not yet 
available; nor are they likely to become available any time soon for some of the more 
complex processes where giant i^'-factors have been observed (e.g. some observables in 
pp — )■ Wbb |13l Alternatively one could simply try to avoid observables like ptji 

and -f^Tjets in inclusive event samples. For example, with additional cuts on the vector- 
boson momentum or a second jet, refs. [121 [U] showed that the i^'-factors are significantly 
reduced. However, given the many analyses that are foreseen at the LHC, it is likely that 
at least a few will end up probing regions where giant i^-f actors are present. 

To understand how else one might address the problem of giant i^-factors, one can 
observe that in our Z+jet example, the bottleneck in obtaining a NNLO prediction is the 
inclusion of the two-loop 2 — )■ Z + 1 parton contributions and proper cancellation of all 
infrared and coUinear divergences. Yet the two-loop (and squared one-loop) contribution 

^This differs from double electroweak (EW) logarithras, which involve terms like aEwln^pt/mz , and 
are usually much smaller. Examples do exist of "giant" EW effects when tagging flavour 211. 

^Part of the enhancement at high pt also comes from the fact that one can have qq — )■ qq scattering 
that emits a Z, whereas the qq partonic channel does not contribute at LO. 
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will have the topology of diagram A in fig. |2] and should not be responsible for the 
dominant part of the NNLO correction, which will instead come from diagrams with the 
topology of B and C, with either an extra QCD emission or a loop. So if one includes tree- 
level 2 — 7- Z + 3 and 1-loop 2 — )■ Z + 2 diagrams (i.e. Z + 2 jets at NLO) and supplements 
them with even a crude approximation to the two-loop 2 — )■ Z + 1 result, one that suffices 
merely to cancel all divergences, then one should have a good approximation to the full 
NNLO result (a related observation has been exploited to obtain approximate NNLO 
results for high-pt J/ip production in [23l|). 

The purpose of this article is to develop a general method for obtaining such rough 
estimates of missing loop corrections. Our approach, called LoopSim, will be based on 
unitarity. After explaining how it works in section [21 and outlining a secondary "reference- 
observable" approach for control purposes, we will test the method by comparing its 
results to full NNLO predictions for lepton-pf spectra in Drell-Yan production in section HI 
apply it to our Z+jet observables in section [5] and finally, in section [6], examine whether 
it can be of use even in the absence of giant i^-factors, specifically for a number of dijet 
observables. 

2 The LoopSim method 

The main ingredient of the LoopSim method is a procedure for taking a tree-level event 
with n final state particles and supplementing it with a series of events with n — 1 particles 
(approximate 1-loop events), n — 2 particles (approximate 2-loop events), etc., such that 
the sum of the weights of the full set of events is zero. This "unitarity" property will 
ensure that all the soft and collinear divergences of the tree-level matrix elements will 
cancel against identical divergences in the simulated loop contributions. 

An outline of the procedure is given in fig. [21 Given a tree- level input event (a), the 
first step is to interpret it as a sequence of emissions (as if it had been produced by a 
parton shower), so that for example (diagram b) one can view particle 2 as having been 
emitted from particle 1, and particle 4 as emitted from the beam. The attribution of an 
emission sequence can be performed with the help of a suitable sequential-recombination 
jet algorithm and will be most meaningful in the limit that emissions are strongly ordered 
in angle and energy. The next stage is to decide which particles refiect the underlying hard 
structure of the event. If the event structure at the lowest possible order is that of a 2 — )■ 2 
scattering, then one should identify two outgoing "Born" particles. The Born particles 
will remain present in all the approximate "loop" events that are generated. They are 
represented as thick red lines in diagram (c). Again this step is most meaningful when 
all non-Born emissions are soft and collinear. 

One then generates a set of simulated "1-loop" events by finding all ways of recom- 
bining one emitted particle with its emitter, diagrams (d,e). Each such "1-loop" event 
comes with a relative weight of —1 compared to the tree-level diagram. Similarly the set 
of simulated "2-loop" events is obtained by finding all ways of recombining two emitted 
particles with their emitter (s) (diagram f), each with relative weight +1; and so forth 
down to events where only Born particles remain (in fig. [3] this is already reached at the 
two-loop level). Note that the loop-diagrams drawn in fig. [3] are not intended to represent 
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(e) 2nd output 1-loop event (f) Output 2-loop event 



Figure 3: Sketch of the LoopSim procedure as apphed to a tree-level event (a) with 4 
outgoing particles (numbered) and the beam (horizontal line); diagram (b) shows the 
attribution of the emission sequence, (c) the identification of the Born particles (thick 
red lines), and (d)-(f) the resulting "looped" diagrams. These diagrams are relevant in 
approximating next-to-next-to-leading corrections to a process whose LO contribution 
has a 2 — )• 2 structure. 



the actual Feynman diagrams that would be relevant at 1 and 2-loop level. Instead they 
indicate the way in which we have approximated the loop divergences, as the unitarising 
counterparts of the divergences that appear for each emission in the soft and collinear 
limits. 

Given the above procedure for unitarising tree-level events, we shall see that it is then 
straightforward to extend it to event sets that also include exact loop diagrams. 

2.1 The tree- level pure glue case 

We start by examining the LoopSim procedure in the simple case of purely gluonic tree- 
level events. This will suffice to introduce most of the relevant concepts. Section 12.21 will 
then discuss some of the additional issues that arise for events with quarks and vector 
bosons, while the handling of events sets that include exact loop diagrams will be left to 
section 12.31 

It is helpful to introduce some notation: Firstly, b is the number of final-state particles 
present in the lowest relevant order (i.e. the number of final-state "Born" particles). For 
instance 6 = 2 if considering higher-order corrections to dijet events, as in fig. [31 
represents a generic event with n final state particles. So the starting event of fig. E] 
would be labelled E4. Finally, f/f will be an operator that acts on an event En and 
returns all the events at / loops obtained from En using the LoopSim method. For 
instance, fig. |3li,e represents the action of Ui^^ on the input i?4 event (a). 

The central part of the LoopSim method involves the construction of the operator 
Ui acting on En for all / = . . . n — 6 (/ < n — 6 because the number of real final state 
particles cannot be smaller than that of the lowest order event). 
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2.1.1 Attribution of structure to events 

Recall that the primary function of the LoopSim method is to cancel the divergences 
that appear in the soft and collinear limits. In these limits, events can be interpreted 
as stemming from a sequence of probabilistic (parton-shower) type 1-^2 splittings of 
some original hard Born particles. The knowledge of the splitting structure will help us 
generate loop events to cancel the divergences. 

The attribution of a branching sequence is most easily performed using a sequential 
recombination jet algorithm (and is inspired by the CKKW matching procedure [21]). We 
will use the Cambridge/ Aachen (C/A) algorithm [23122], which has important advantages 
over the kt algorithm when dealing with nested collinear divergences (avoiding "junk" jets 

As a first step, to each of the i = 1 . . .n particles in the event En, we assign a unique 
"identity" index li = i. 

We then run the Cambridge/Aachen (C/A) algorithm on the event. It repeatedly 
clusters the pair of particles that are closest in angle, i.e. with smallest dij = AR'fj/Rl^ 
where ARfj = {yi — yjY + (0j — 0j)^ is the usual squared angular distance in the {y, (j)) 
plane, and i?LS is a free parameter, the LoopSim radius. The C/A algorithm continues 
until all the dij > 1, at which point the remaining particles are deemed to cluster with 
the beam. 

An ij — !■ k clustering in the C/A algorithm can be reinterpreted as a — )■ ij splitting. 
The C/A algorithm does not distinguish in any way between i and j. However, in the 
soft limit, say ptj <^ pu, rather than viewing k as splitting to i and j, it is a better 
reflection of the divergent structure of the amplitude to view k as having emitted a soft 
gluon j. Then i is nothing other than particle k with some small fraction of its energy 
removed. To account for this, in an ij — )■ k clustering, if ptj < pu, then we declare that 
the "identity" Ik of particle k should be the same as that of particle i, Ik = li- Also we 
record li as being a "secondary emitter" and remember that the object with identity Ij 
has been emitted from the object with identity /«. (Exchange i -H- j if pu < ptj)- This 
is represented in flg. [3]d by the fact that particle 1 is a straight line, off which particle 
2 has been emitted; the identity of the 1 + 2 combination is J1+2 = Ji = 1. For an iB 
clustering, we record /j as having been emitted from the beam. 

The next step in attributing structure to the event is to decide which event particles 
should be viewed as the Born particles, i.e. which particles are responsible for the hard 
structure in the event. Inspired by the original formulation of the Cambridge algorithm 
^25j , for every ij — )■ k recombination we assign a fcf algorithm type hardness measure hij = 
mm{pf^,pfj)ARjj/Rlg [221 [3Oj0 For every beam recombination, we assign a hardness 
hiB = Pti- 

We then work through the recombinations in order of decreasing hardness. For an 
ij — )■ k recombination (or k — )■ ij splitting), assuming i is harder than j, we mark Ik = li 
as a Born particle. If fewer than b particles have already been marked as Born particle, 
we also mark Ij as a Born particle. For an iB recombination, we mark /j as a Born 

•^AU jet clustering in this article is carried out using Fast Jet p71 128| . 

"'In the results shown later, we actually used hij = min{p^^,Pfj)ARfj, which however is identical for 
our default choice of i?LS = 1- 
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particle. This is repeated until b particles have been marked as Born (a particle may be 
marked more than once; in such a case its marking counts only once). As an example, in 
fig. [31 the hardest recombination will be between particle 3 and the beam, so particle 3 is 
marked as a Born particle. The next hardest recombination will that of (1 + 2) with the 
beam. Therefore we mark /(1+2) = /i = 1 as a Born particle. This exhausts the number 
{b = 2) of Born particles that need to be marked. 

At the end of the above procedure, every particle will have been marked as emitted 
either from the beam or from another particle, and some particles will also have been 
marked as secondary emitters and/or Born particles. Thus in figure |3l particle 1 is labelled 
as having been emitted from the beam, it is a secondary emitter and a Born particle; 
particle 2 is labelled as having been emitted from particle 1; particle 3 is a Born particle, 
emitted from the beam; and particle 4 was emitted from the beam. The structure that 
we attribute is of course physically unambiguous only in the presence of strong ordering 
of emission angles and energies. However, as we shall argue in section 12.41 the mistakes 
that we make for non-ordered configurations should have a small impact for observables 
with giant i^'-factors. 

2.1.2 Constructing virtual (loop) events 

Once every particle is labelled in an event En, one can compute the result of Ui{En), 
which is a set of events En-i- For an event En with respectively b Born particles and 
non-Born secondary emitters, we define 

v = n-{b + ns), (1) 

to be the maximum number of particles that will be allowed to become virtual in a given 
event. It is obvious that Born particles will not become virtual. Additionally, secondary 
emitters will also not become virtual. To understand why, consider the event 





1 




2 



(2) 



in which particle 3 is a secondary emitter, since it emitted 4. There is a divergence for 
4 to be coUinear to 3 only if 3 is a final-state particle. If instead 3 is made virtual, then 
the divergence for emitting 4 no longer exists (there is no divergence for emission from 
internal lines in a diagram), so that the weight for the diagram in which 3 is virtual would 
be the weight of the tree- level diagram times a small coefficient e <^ 1. This simplest 
way of accounting for this is to approximate e = and thus not generate events in which 
secondary emitters are made virtual (a more detailed discussion is given in appendix O). 

Having understood which particles can be made virtual, the operator Ui, when applied 
on an event En, generates all the (^) diagrams in which / particles become virtual. For 
the virtual events to cancel the infrared and collinear divergences that appear in the tree- 
level diagram, we need an infrared and collinear (IRC) safe procedure to make particles 
virtual. For instance, the divergent weight of an event with two collinear partons i and j 
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has to be cancelled by that of corresponding virtual event {j makes a loop over i) when 
computing the distribution of any IRC safe observable; and two coUinear partons, if not 
virtualised, have to remain coUinear when another particle becomes virtual. 
There are two ways for a particle j to make a loop: 

• If it is labelled as clustering with particle i, then one has to spread the momentum 
of particle j over i and all the particles that are labelled as clustering with it but 
which were emitted after j according to the C/A clustering sequence (i.e. at smaller 
angle). The exact procedure is explained in detail in appendix |X1 and is designed to 
ensure that the recombination maintains any collinearity properties of non-looped 
particles and is invariant under longitudinal boosts. When j is the only particle 
that clusters with i, then the procedure becomes equivalent to adding the momenta 
of particles i and j , Pk = Pi + Pj, and then rescaling the momentum p^ such that 
its mass is set to 0, while leaving its transverse components Px, Py and its rapidity 
unchanged. 

• If particle j is labelled as clustering with the beam, then when it is "looped" it is 
simply removed from the event. Note that looping particles with the beam is less 
trivial than it may seem at first sight, because of an interplay with factorisation 
and the PDFs. Nevertheless it can be shown, appendix [B], that for particle types 
that are included in the PDFs it does make sense to loop them. A pt imbalance 
will result from the looping of particles with the beam, and so after all loops have 
been made, we apply a transverse boosts to all remaining event particles, conserving 
their rapidities, so as to bring the total transverse momentum to zero (again, see 
appendix El). 

There is some arbitrariness to our procedures for producing physical kinematics in the 
looped events. One avenue for future work would be to examine the impact of making 
different choices. 

The operator Ui has the following properties 

Ul = \, U^{En) = if / > . (3) 

If Wn. is the weight of event En, then each of the events generated by the Ui{En) operator 
has a weight 

Wn-l = (-1)'W„ . (4) 

Once all the Ui{En) have been calculated for I = . . .n — b, one has to combine them in 
order to subtract all the soft and collinear divergences that appear in the calculation of 
En and the virtual diagrams generated from it. This is done by the operator f/y, which 
is defined as 

V 

u^^Y^^r (5) 

1=0 

It generates all the necessary looped configurations that have the same order in as 
the original tree-level diagram. It is straightforward to see that the total weight of the 
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diagrams obtained from the operator is 0. Indeed, if we apply it to an event En whose 
maximum number of virtual particles is v, we get 

for t> > 0. 

We note that the above procedure for approximating loop diagrams does not generate 
the finite terms needed to cancel the scale-dependence of lower-order diagrams. While it 
would be straightforward to include such terms, we believe that in the absence of full loop 
calculations, not including them helps ensure that the standard procedure of variation of 
renormalisation and factorisation scales is more likely to provide some form of reasonable 
estimate of the uncertainties on our results. 

2.1.3 Some examples 

In order to illustrate the action of the operator Uj', we give below some simple examples 
in the pure glue case. In each of these examples, only the Born particles are labelled with 
numbers 
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6=2 
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(7f) 



Eq. f l7ap gives an example of singly-looped configurations ("1-loop diagrams") generated 
by LoopSim when studying the 2 — )■ 4 contributions to QCD dijet production. Eq. flTbj) 
shows the "2-loop diagrams" generated from the same event. The next equation shows 
what happens if we add one more particle to the final state. If, eq. (l7dl) . we now set the 
number of Born particles for the same event to be 3, we obtain only one 2-loop diagram 
instead of three, as represented in eq. (l7d|) FI Finally, the last two examples of eq. ([7]) give 
a case with a splitting: the emitter is not looped, even if it is not a Born particle. 
We also give a few examples of the action of the ^7^ operator: 
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In the last case, only one particle can become virtual because there are two secondary 
emitters which cannot be looped. 



2.2 Treatment of flavour within LoopSim 

Let us now examine some of the issues that arise if we are to extend the LoopSim method 
to processes with quarks and vector bosons. 

We start with quarks and consider the situation depicted in fig. HI In this case, 
applying the C/A algorithm as in the previous section will lead to the recombination of 
the two quarks qi and q2, which is clearly not physical. If flavour information is available 

^One might reasonably be surprised by this: after all, the result for the exact two-loop diagrams is 
independent of any choice of number of Born particles. The point is that if one studies soft and coUinear 
corrections to the 3-jet cross section, then for the events in eq. ([7c)) where the particle labelled 3 in 
eq. (I7dl) is virtual, the event will resemble a two-jet event and so not pass the 3-jet cuts. However if 
one studies the 3-jet cross section in a kinematic region where the cuts allow one of the jets to be much 
softer than the others, then to obtain sensible results it becomes necessary to use 6 = 2 and include all 
diagrams on the right-hand side of eq. ((7c|) . 
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Figure 4: Example of an event where two quarks qi and q2 may get recombined by the 
C/A algorithm. 

for events, then one can veto on such a clustering, for instance by defining the clustering 
distance dqq between two quarks to be infinite. As discussed in [31] such a modification 
alone is not sufficient to systematically guarantee sensible treatment of flavour in jet 
clustering. Refs. [311 132] have both discussed the further modifications needed in the case 
of the kt algorithm. A proper handling of flavour within LoopSim might seek to extend 
those modifications to the C/A algorithm. However, neither of the NLO programs that 
we use, MCFM and NLOJet++, provide information on particle flavours, so we defer 
such modifications to future work and just maintain the dij = ARfj/Rf^^ distance for 
all partons. For observables that are not flavour-sensitive this should not be a major 
drawback, given the observation [31] that divergences associated with the mistreatment 
of flavour are strongly subleading. Were we to be interested in heavy tagged quarks, more 
careful treatment might well be needed. Note that there are also subtleties related to 
flavour and PDFs, discussed in appendix iBl 

What about non-QCD particles, specifically vector bosons? Let us examine the case 
of Z bosons. A Z can be emitted from quarks or antiquarks and we would like this to 
be reflected when establishing the approximate emission sequence, because if the Z has 
been emitted from a quark, then that quark is a secondary emitter and should not be 
looped. In other cases a Z boson may be the hardest isolated object in an event. Then 
it is to be considered a Born particle. On the other hand we won't necessarily wish to 
consider diagrams where a Z boson is looped, because they would represent electroweak 
corrections, not QCD corrections. 

One issue in dealing with electroweak particles is that they are not emitted from 
gluons. If one could distinguish between quarks and gluons, then this could be accounted 
for during the C/A clustering, by defining the distance dgz between a Z and a gluon to be 
infinite. Since we will not know which partons are quarks or gluons, we adapt Frixione's 
isolation procedure p3j to decide if a Z boson relatively close in angle to a parton i is 
likely to have been emitted from i. More precisely, if 



then we define diz = AR'^^/Rl^, otherwise diz = oo. When recombining i and Z into a 
particle k, then the identity index Ik is set equal to li (a quark and a Z give a quark). 
Our procedure means that a Z that is very collinear to a parton is always considered to 
be emitted from that parton — this makes sense because such configurations are much 
more likely to occur when the parton is a quark. In contrast a soft parton in the general 




(9) 
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vicinity of a Z is not clustered with the Z, which is sensible given that most soft partons 
tend to be gluons. Finally, for a recombination between a parton i and Z, we define the 
hardness of the branching hiz a^ 



hiz = min(p|, + ml)—^ , (10) 

while a recombination of a Z with the beam has a hardness 

hzB=Ptz + mz- (11) 

The latter means that for an event with just a parton and a recoiling Z boson, the parton's 
beam hardness will always be lower than the Z's, implying that for 6 = 1 it is the Z-boson 
that will be the single Born particle, as should be the case, at least when the parton has 
Pt -C mz, i.e. in the kinematic regime that dominates the total cross section for the Z. 

Once a structure has been assigned to an event with a Z boson, the next question is 
that of the looping procedure. When looping partons it remains identical to before, with 
just a small extension of the recoil procedure in order to deal with decay products of the 
Z boson (see appendix . In the situations where the Z is not a Born particle (it is 
never an emitter), straightforwardly following the procedure of section [2.1.21 one would 
deduce that one should loop the Z as well: 



y 2 l^z y 2 2^ 
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(straight lines are partons, either quarks or gluons). The rightmost diagram, with the 
looped Z, is not, however, a QCD loop diagram: it is an electroweak loop correction to 
a multijet event. The LoopSim procedure does not aim to reproduce electroweak loop 
corrections (though in this case it might be a reasonable approximation). Furthermore, 
in any analysis that tags on Z bosons, such a diagram would not be tagged and so would 
not contribute. Thus, although the LoopSim procedure naturally generates events with 
looped Z bosons, events like the rightmost diagram of eq. ffT^ are simply to be discarded. 

For events with bosons, the same procedure can be used as for Z's. Note, how- 
ever, that while the "looped" Z events may give a reasonable approximation to actual 
electroweak loop diagrams, looped events will not. This is because W-boson emission 
changes quark flavour: consider a tree- level diagram ud — > bbW^, with the W"*" emitted 
collinearly off the incoming d, converting it into a u. The LoopSim procedure would give a 
"loop" diagram ud — )■ bb, with the W"*" looped. However no such loop diagram exists and 
the correct loop diagram would instead involve uu — )■ This is closely related to the 



^This will be true in a future version of the code, but currently hiz = m.iii(p^^,p^,^)AR^^. 

similar problem would appear to exist with the QCD diagram ug — >■ bbu, with the incoming gluon 
splitting collinearly to give uu, and the outgoing collinear u being looped. Here, however, we are saved 
by the interplay between LoopSim and PDFs, as discussed in appendix |B] It is crucial in this respect 
that all flavours that get looped are included in the PDFs. 
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phenomenon of Bloch-Nordsieck violation [31] that is found when considering electroweak 
double logarithms. Since we in any case discard events in which electroweak bosons are 
looped, this should not be a problem for the practical use of LoopSim in events with 
bosons. 

2.3 Merging NLO calculations and beyond 

Before explaining how we merge exact higher orders calculations, let us mention how we 
use the LoopSim method in practice on tree-level events at several different orders. We 
introduce the notation XtQln^LO to denote an approximation to the N^LO cross section 
for producing X, with all loop terms estimated through the LoopSim procedure. It is 
obtained by applying the operator to all tree-level diagrams that can contribute up 
to N^LO. For instance, one can write 



Notice that f/^(Z@LO) = Z@LO and ?7^(Z+j@L0) = Z+j@LO. The terms ?7^(Z+j@L0) 
and ?7y (Z-|-2j@L0) simulate up to one- loop corrections, and [/y(Z+3j@L0) simulates up 
to two-loop corrections. 

Now let us see how things work beyond tree-level accuracy. We define En,i to be 
a generic event at I loops (exactly calculated) with n particles in the final state. We 
first consider the case where only one-loop corrections are computed exactly, so that we 
have tree- level events Enfl and exact one- loop events En-i^i- As before we can apply the 
unitarisation operator to the tree-level events, U^{En,o)- However, since we now include 
exact 1-loop contributions, En-i,i, we must remove the approximate 1-loop contributions 
Ui{Enfi) that are contained in U^{Enfi)- This alone is not sufficient, because among the 
extra contributions from the exact 1-loop terms, there will be pieces that are finite for a 
given (n — l)-parton configuration, but that can lead to divergences when integrated over 
the (n — l)-parton phase space. To cancel these extra divergences, we should introduce 
additional approximate higher-loop contributions, which can be obtained by applying the 
unitarisation operator f/y to the difference between the exact and approximate one-loop 
terms. So, rather than including just events -En-i,! and subtracting Ui{Enfl), we include 
events U^{En-i^i) and subtract (U^{Enfl)) ■ It is convenient to express this through a 
new operator ^ such that 



exact 1-loop result is to be included. The action of f/y ^ depends on the number of loops 
already included in the event on which it operates: we subtract the one-loop contribution 
returned by LoopSim only in tree-level events. With this notation, one can compute the 



Z+j@nLO 
Z-|-j@nnLO 



Z@nLO 



?7^(Z@L0) + [/v^(Z+j@LO) , 

;7^(Z+j@L0) + f/^(Z+2j@L0) , 

[/^(Z+j@LO) + f/^(Z+2j(aL0) + f/^(Z+3j@L0) . 



(13a) 
(13b) 
(13c) 
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higher order corrections to eqs. (1131) to one-loop accuracy, 

Z@nNLO = Z@NLO + f/v^_i(Z+j@NLO„„iy) , (15a) 
Z+j@nNLO = Z+j@NLO + f/^^i(Z+2j@NL0„^iy) , (15b) 
Z+j@nnNLO = Z+j@NLO + f/^ ,(Z+2j@NL0„,iy) + f/^ i(Z+3j@NL0„,iy) , (15c) 

where the "only" subscript on Z-|-nj@NLOoniy means that we take the highest order that 
contributes, i.e. here a^~^^aEw, since the LO, a'^a^^w, piece of Z-|-nj@NLO, is already 
taken into account in the Z+{n — l)j@NLO contribution. This implies that one should 
use consistent renormalisation and factorisation scale choices across all different orders of 
the calculation. Note that in eq. (|T5|) we have introduced the notation n^N'^LO to denote 
an approximation to the N^+'^LO result in which the p highest loop contributions have 
been approximated with LoopSim. 

The extension of the procedure beyond one-loop accuracy is simple. For instance, at 
two-loop accuracy, one has to subtract the approximated two- loop contribution [/^(-^n.o)" 
Ui (Ui{En,o)) in eq. fll4al) . and the other approximated two-loop contribution Ui{En,i) 
in eq. fll4bl) . giving 

U^^,{E^,o) = U^{E^,o) - (U'AEnfi)) - f/v KiEnfi) - f/f {UKEnfi))] , (16a) 

[/^,2(^n-l,l) = f/^(K-l,l) - (^7i^(S„-l,l)) , (16b) 
f/v,2(^n-2,2) = f/v'(^-2,2). (16c) 

Therefore, once Z-|-j@NNLO is calculated, one may compute for instance 

ZOnNNLO = Z@NNLO + f/v^_2(Z+j@NNL0„„iy) . (17) 
To be complete, let us mention the generalisation of our procedure to m-loop accuracy 

m—l 

c/^,j£;„_,.,) = f/v'(i?„-M) + E(-i)' E ° K ° ■ ■ ■ ° ul{E^-i,i) . (18) 

i=i h,---,h > 1 

h + . . . + Ij < m — I 

We noted at the end of section 12.1.21 that the plain LoopSim procedure does not 
generate the finite terms needed to cancel residual scale dependence from lower orders. 
With the introduction of the exact loop contributions, those finite terms do now get 
included. Thus for a given number of exact plus simulated loops, as we increase the 
number of exact loops, we should expect to see reductions in scale uncertainties. 

2.4 Expected precision of the method 

Let us briefly explain why the LoopSim method is expected to work in the presence of 
giant i^'-factors. We consider an observable A computed respectively at NLO and nLO. 
We define -R'nlo ^^"^^ ^^^^ 

and we assume that 3> 1. This huge ii'-f actor may come from logarithmic en- 

hancements in the real NLO diagram or the appearance of new scattering channels in the 



15 



perturbative expansion. The computation of ct^q gives the exact real part of the NLO 
calculation as well as the divergent terms of the virtual correction. Therefore 

4?o - = O , (20) 

where, in writing OlagC^), we mean that the term missing in the nLO calculation, the 
finite part of the 1-loop correction, is not especially enhanced. This leads to 

4?o = 4fo [l + O I I . (21) 

The relative difference between the approximate and exact NLO calculations is thus 
suppressed by the inverse K-factor. 

Next, consider nNLO accuracy. The difference between 4nlo "^nnlo comes from 
the parts of the two-loop corrections that are finite and associated with the LO topology, 
so that they should be free of the enhancements that led to the large NLO i^'-factor. This 
implies 

41o - 41o = O • (22) 

If we define K^-^q such that ct^^^q = -^^nnlo'^lo ' then we can write 

41o = 41o + I I . (23) 



-"^NNLO. 



If -R'nlo ^ 1' Cine can expect -R'nnlo ^ ^ ^^lo. 



3 The reference-observable method 

Given the novelty of the LoopSim method, it is useful to have an alternative way of 
estimating the size of the NNLO contributions that we will approximate with LoopSim. 
Here we outline such an alternative method, which, though less flexible than the LoopSim 
approach, will provide a valuable cross-check and help us build our confidence in results 
of the LoopSim method. 

Let us explain it for observables in the Z+j process. Our aim is to estimate ct^^^q for 
some observable We assume that we have a reference observable which is identical to 
the observable A at LO. For instance, one might consider ref = pt^z and A = ptj. We can 
write the NNLO Z+i prediction for A in terms of the NNLO prediction for the reference 
observable plus the NLO Z+2j difference between A and the reference cross section 

4+j<aNNLO = 4+j@NNLO + i'^'^^^ " o-^''''''^)z+j@NNLO , (24a) 
= 4+?@NNLO + (^^^^ - cr^'''%+mNLO . (24b) 



^More precisely, the cross section for the observable A to pass some given cuts; it is only for 

brevity that we use here the somewhat inaccurate shorthand "cross section for observable A" . 
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The second equality is possible because 2-loop NNLO corrections to Z+j have the topology 
of Z+j at LO. Therefore, their contributions to the observables A and ref are identical 
and cancel in the difference in eq. fl24ap . 

If we have reason to believe that the perturbative expansion for the reference ob- 
servable converges well, we can conclude that CTz+jeNNLO ~ ^z+j@NLO genuinely a small 
correction. Then 

(A) (ref) / (A) (rcf)\ fnr:\ 

'^Z+j@NNLO-^Z+j@NLO + l'^ 'jz+2j@NLO , l^Oj 

i.e. we approximate the NNLO distribution for A in terms of the NLO distribution for the 
ref observable and a NLO calculation for difference between the A and ref distributions, 
both of which are exactly calculable. The missing part is suppressed by a relative factor 
1/K^'^\ as for the LoopSim method. For Z+j, one can see from fig. [T]that pt^z seems to 
be an acceptable reference observable for ptj and -ffrjets- 

In the sections that follow we shall, for brevity, refer to the RHS of eq. (1251) as "ref. 
TT-NLO" even though it does not quite adhere to our the meaning of riNLO as set out in 
section 2, i.e. in terms of the specific sets of tree-level and loop diagrams that are included 
exactly. 

4 Validation: comparison to DY at NNLO 

The cross section for the Drell-Yan process is known with exclusive final states up to 
NNLO accuracy [551 EH]- Above a certain value of lepton transverse momentum, one 
finds giant corrections to the lepton pt spectra when going from LO to NLO and large 
ones from NLO to NNLO. This gives us an opportunity to directly test the performance 
of the LoopSim method by comparing its nNLO results to exact NNLO spectra for lepton 
pair production. 

Before examining nNLO results, it is useful to compare nLO with NLO. If they are 
in reasonable agreement for some observable, then that serves as a first indication that 
the LoopSim estimate of missing loop corrections is sensible for that observable. 

Fig. Ogives the comparison of the nLO, NLO and LO results for the production of 
an e"^e~ pair within the mass window of 66 < rrie+e- < 116 GeV at a proton-proton 
centre of mass energy of 14 TeV. The left-hand plot shows the cross section differential 
in the transverse momentum of the harder of the two leptons. The right-hand plot 
gives the corresponding K factor with respect to LO. The results were obtained with 
MCFM 5.3 [an 119], with its default set of electroweak parameters and NNLO MSTW2008 
parton distribution functions. The uncertainty bands in Fig. [5] correspond to varying the 
renormalisation and factorisation scales /Ir = /ip by a factor of | and 2 around a default 
choice of mz- In the nLO result we fixed the value of the LoopSim radius parameter 
to be i?LS = I5 which naturally places interparticle and particle-beam clustering on the 
same footing (though the nLO result here is actually independent of i^LS, because there 
is at most one isolated QCD parton in the final state). 

There are three relevant regions of transverse momentum in fig. |5l For pt,max ^ 
(low Pt) the distribution is dominated by on-shell Z-bosons and its shape is governed by 
the angular distribution of the Z decays in their centre-of-mass frame. The peak close 
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Figure 5: Comparison between nLO results from LoopSim (with i?Lg = 1) and exact NLO 
results for the Drell-Yan process. The left-hand plot shows the transverse momentum 
spectrum of the harder lepton, while the right-hand plot gives the corresponding K factors 
w.r.t. LO. The uncertainty bands were obtained by varying /Xr = yUp by a factor of | and 
2 around a default choice of mz- 



to |mz corresponds to Z-bosons that decay in a plane at right-angles to the beam. For 
l^z ^ Pi,max < 58 GeV (intermediate pt), the LO distribution comes from Z-bosons 
that are off shell, which allows the pt of the lepton to be larger than |mz. The narrow 
width of the Z causes the distribution to fall very steeply. The 58 GeV upper edge of this 
region is a consequence of our cut on rrie+e- < 116 GeV. Above 58 GeV (high pt) the LO 
distribution is zero. 

In the low pt region, the NLO correction is moderate and negative. There is no strong 
reason to believe that the LoopSim method should work here, but it turns out that 
the nLO result reproduces the structure of the correction, even if its scale dependence 
remains much larger than that of the NLO result (this is because the LoopSim procedure 
does not include the finite terms that would partially cancel the LO scale dependence). 
In the intermediate pt region, we see a "giant" NLO i^-factor. It comes about because 
initial-state radiation can give a boost to the Z-boson, causing one of the leptons to 
shift to higher pt (it becomes the "max" lepton). The spectrum of QCD radiation falls 
much less steeply than the Z-boson lineshape, so this NLO correction dominates over 
the LO result. In this region the exact loop correction, proportional to the LO result, 
becomes almost irrelevant and we see near perfect agreement between nLO and NLO. 
In the high-pt region only the real emission diagrams of Z(QlNLO contribute and nLO 
becomes identical to NLO (both correspond to the Z+j@LO result). Similar results hold 
for the pt^e± distribution, while the pt,min lacks the giant i^'-factor in the intermediate 
region. 

A similar comparison between nNLO and NNLO spectra is shown in fig. Ei The NNLO 
results were obtained with DYNNLO 1.0 OB EHl EH], used with a set of electroweak 
parameters compatible with that of MCFM|j 

^ In its 0{as) and Oi^a"^) contributions, DYNNLO includes among its parameters a cut on the pt of 



18 




Figure 6: Comparison between nNLO results from LoopSim+MCFM (with i?Ls = 1) and 
full NNLO results for the Drell-Yan process from DYNNLO. The left-hand plot shows 
the transverse momentum spectrum of the harder lepton, while the right-hand plot gives 
the corresponding K factors w.r.t. NLO. The uncertainty bands come from varying the 
factorisation and renormalisation scales by factors 1/2 and 2. In the right-hand plot we 
also show the (thin) band related to changing the nNLO /2ls parameter from 0.5 to 1.5, 
at fixed /Ir = /Xp = mz- 



In the low-pi region we find quite good agreement between the nNLO and NNLO 
results (with somewhat larger uncertainty bands for nNLO). Such a result was not 
guaranteed a priori, even if it is not entirely surprising given the reasonable agreement 
that we saw between nLO and NLO. In the intermediate pt region, where the NNLO/NLO 
corrections are substantial, the agreement is excellent. This was expected. At high pt 
the agreement should be exact, and does seem to be, within statistical fluctuations. The 
dependence on i?LS (shown in the right-hand plot) has been estimated by varying its 
value from 0.5 to 1.5. The effects are small. 

Finally, we note that similar features and a similar level of agreement between nNLO 
and NNLO are to be found in the pt,mm and pt^e± distributions. 



5 Results for the Z+jet process 

In the previous section, we studied the Z production process and showed that our pro- 
cedure correctly reproduces the pt distribution of the hardest lepton at NNLO, even, 
unexpectedly, in regions where the ii'-factor is not large. In this section we study the 
Z+i process, whose NNLO cross-section is not known yet, but which leads to giant K- 

the Z boson. The cut is applied to both real and virtual terms and its impact should vanish as it is taken 
towards zero. It is, however, required to be non-zero for the numerical stability of the MCFM Z+j NLO 
calculation that is among the components of DYNNLO. We set the cut equal to 0.1 GeV in the 0{as) 
term and to 1 GeV in the 0[af) term. A related 1 GeV cut was placed on the 0(^a^) piece of the nNLO 
result (while none was used at 0{as))- The impact of the 1 GeV cut is small but not entirely negligible 
close to the peak (where, physically, NNLO should in any case be supplemented with a resummation) . 
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factors at NLO for some observables as explained in the introduction. Therefore, their 
NNLO contributions are expected to be accurately described by the LoopSim method. 
Throughout this section we use MCFM 5.7, including the Z+2j process at NLO [40j, with 
the NLO CTEQ6M PDFs. We will take three different values for the renormalisation and 
factorisation scales: /i^ = /Xp = |/xo, /xq and 2/io, with 

/^o = -^Z"^! (26) 

where ptji is the transverse momentum of the hardest jet. At high pt, this scale choice 
should be quite similar to that used in [41] and has the same pt scaling as those in 
[m [15] . The i?LS uncertainty is measured at /Xr = fip = fio using three different values 
for it: i?LS = 0.5, 1, 1.5. 

In addition to the 3 observables shown in the introduction, z, ptji and -f^Tjets = 
Yl'j^iPtjiy we will also consider 

-f^T,tot = -f^Tjets + Pt,z ■ (27) 

We only include events for which ptji > 200 GeV. 



5.1 Validation at nLO 

As a first investigation of the performance of the LoopSim method, let us examine how 
the nLO approximation compares to the full NLO result. Fig. [7] shows the i^'-factors for 
the fiLO and NLO predictions, with uncertainty bands from scale and -Rls variations. 

In the upper- left plot, one sees that the nLO prediction for the pt^z distribution gives 
a somewhat smaller /^-factor than the NLO result. We interpret this as being because 
certain genuine loop effects are not taken into account by the LoopSim method, for 
example those related to threshold logarithms, which depend crucially on the factorisation 
scheme of the parton distribution functions. The nLO result does, however, reproduce 
the Pt dependence of the i^-factor, i.e. the dip towards pt = 200 GeV. This dip arises 
because of the requirement in our event selection that there should be at least one jet with 
Pt > 200 GeV. At LO this induces a step-function in the pt^z distribution at 200 GeV. At 
NLO, soft and collinear emissions smoothen out that threshold and the nLO calculation 
correctly reproduces the resulting interplay between real and virtual terms. 

In the three remaining plots of fig. U\ for Ptji, -f^Tjets and -ffr.tot, all of which have giant 
i^'-factors, one sees good agreement between the nLO and NLO results. This is because 
the dominant NLO contribution comes from events in the B and C-type configurations 
of fig. m for which there is no corresponding QCD loop correction. The LoopSim method 
merely serves to cancel the divergences that arise from soft and collinear emissions off 
A-type configurations and these are not dominant overall. 

The i?Ls dependence, also shown on these four plots, only comes from 1-loop events 
generated by LoopSim. Therefore, for an observable A studied in Z-|-j@nLO with two 
different values Rq and Ri for i?Lg, one can write: 

AA) _ (A) _ (pt,z) _ (Pt,z) /r,oN 

°Z+j@nLO,Ri °Z+j@nLO,iio ~ "^Z+jOnLO,/?! Z+j@nLO,iio \'^°) 



20 





Figure 7: Comparison of the nLO/LO fC-factor with the NLO/LO fT-factor, together 
with their scale and i?LS uncertainties for four observables in the Z+jet process. 
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Figure 8: Comparison of the riNLO/LO fC-factor with the NLO/LO i^'-factor, together 
with their scale and i?LS uncertainties for four observables in the Z+jet process 



as long as A coincides with pt^z at LO (it does for each of Ptji, -f^Tjets and |i?T,tot)- This 
means that the absolute uncertainty due to i?LS is the same for A and pt^z- Therefore, 
the relative uncertainty due to i?LS is expected to be roughly inversely proportional to 
the i^-factor for A, in analogy with the discussion of sec. 12. 4[ This explains why the i?LS 
dependence (solid cyan band) looks significantly smaller for Ptji, -f^Tjets and -ffr,tot than 
it does for pt^z plot. 



5.2 Results at nNLO 

Results at nNLO are given in fig. [HI In the case of pt^z the result is similar to the NLO 
result, and the scale uncertainties remain largely unchanged. In other words, since Z+2j 
topologies do not dominate the high-pf z distribution, adding NLO corrections to them 
(i.e. nNLO Z+j) makes no difference either to the result or to the uncertainties. We have 
also shown the dependence on the choice of R in the LoopSim procedure. It is smaller 
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Figure 9: Comparison between the approximate NNLO/LO i^'-factor calculated using 
respectively the LoopSim and the "reference-observable" method for ptji and if-rjets- As 
a reference observable we have used the differential cross section for pt^z- 



than the scale dependence. 

The ptji distribution gets a correction that is just within the NLO uncertainty band, 
with nNLO uncertainties that are about half the size of the NLO band. Adding in the 
77.NLO term has made a real difference. This is precisely what we expect: the observable 
is dominated by Z-|-2-parton configurations, and these were only present at tree-level in 
the NLO Z+i calculation. Our use of fiNLO provides the additional 1-loop Z+2-parton 
and tree-level Z-(-3-parton configurations that come with NLO Z+2j accuracy. 

Given the improvement in scale uncertainty, we need to ask whether the uncertainty 
due to -Rls variation might somehow eliminate part of this benefit. It is, however, small. 
The reasons are similar to those given around eq. (l28i) . 

The -f^Tjets and Hrp^tox. distributions get significant nNLO corrections, with nNLO/NLO 
X-factors of about 1.7 — 2. Absolute scale uncertainties increase slightly compared to 
NLO, but because of the large /^-factor, relative scale uncertainties diminish. At first 
sight, it is somewhat disturbing that the nNLO and NLO uncertainty bands don't over- 
lap. Given the novelty of the LoopSim method, one should therefore ask whether this is 
reasonable and whether there is any way of cross-checking the result. 

A first observation is that since nNLO Z+j is really NLO of the dominant Z+2j 
component, the large nNLO corrections that we see are comparable to an 0{2) K-factor 
for going from LO to NLO in the Z-|-2j prediction. There are many contexts where NLO 
and LO results are not compatible within scale uncertainties, and so it is not unreasonable 
that the same should be seen here. 

Still, we would like to have some more quantitative cross checks that our results 
are sensible. One option is to consider the alternative "reference-observable" method 
presented in section [3], which only makes use of standard NLO calculations to compute 
the approximate NNLO corrections. The comparison between the two methods is shown 
in fig. [9] for -f^Tjets and Pt,ji, where we have taken pt^z as the reference observable. One 
notices near perfect agreement for -ffTjets and very good agreement for Ptji- This gives us 
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some degree of confidence that tfie nNLO LoopSim results provide an accurate description 
of the NNLO behaviour for these observables. 

A second option for cross-checking the large nNLO effects for -ffrjets and -ffr.tot, is to 
examine whether Ht type observables might generally be "difficult". To do so we look 
at them in the case of QCD jet events. 

6 QCD jet events as a testing ground 

We have seen that the nNLO /^-factors for the two effective-mass variables, -f^T.tot and 
-f^Tjets! in Z-l-jet(s) events are about a factor of two above the NLO /^-factor. 

Since NLO is the first order at which we see the dominant "dijet" topology for the 
Ht variables in Z+jet(s), fig. [2]B,C, it might be instructive to establish a correspondence 
with a simpler process, QCD dijet production. Having a NLO Z+j prediction is analogous 
to a LO dijet prediction; and the nNLO Z-|-j predictions should be analogous to NLO 
dijet predictions. NLO cross sections for dijet observables can be calculated exactly 
and therefore we can check whether NLO i^T-f actors of order 2 appear for effective- mass 
observables in pure QCD events. 

We will consider several effective-mass observables: an Ht^u variable, which sums over 
the n hardest jets above some threshold (pf,min = 40 GeV; such a cut is often imposed 
experimentally^"!) 

n 

HT,n= (29) 

iejets with pt,ji > Pt .mill 

where pt^i is the transverse momentum of the i^^ hardest jet. Upper limits on the number 
of jets included in the effective mass are common in SUSY searches We also define 

an effective mass for all jets above the pt,mm threshold, 

Ht = Ht,^ , (30) 

which is similar to the -f^Tjets and i^T, tot observables of section [51 Finally, for completeness 
we will consider the distributions of Ptji, Pt,j2 and the inclusive jet spectrum. All our 
results in this section will be for a centre-of-mass energy of 7TeV, to allow comparison 
to results in the current run of the LHC. 

At LO, the distributions of l-^T.n > 2), ^Ht, Pt,ji, and ptj2 will all be identical. 
The inclusive jet spectrum will have a distribution that is twice as large (because each of 
the two jets contributes). Note that we do not impose any rapidity acceptance limits on 
the jets: though such a cut would have been trivial to include in the LoopSim procedure, it 
would have complicated somewhat the reference-observable approach that we will consider 
at the end of the section. LoopSim results with a rapidity cuts of \y\ < 2 on the jets are 
available from the authors on request. 

Figure [Tor left) shows the distributions for two observables, ^Ht and pt,2 at LO (where 
they are identical) and at NLO, as determined using NLOJet++ ^ gS] with CTEQ6M 
PDFs. A first comment is that receives a NLO i^-factor of order 2, just like the nNLO 

^"^In section [5] we did not apply this kind of cut on the Ht variables; one purpose in applying it here 
is to ascertain whether the large higher-order effects persist even with it. 
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Figure 10: Left: differential cross sections for the •pt^j2 and \liT observables, at LO, 
where they are identical, and at NLO where they have substantially different iiT-factors. 
Right: the NLO i^-factors for the 400 < V j GeV < 500 bin for each choice of variable 
V among the following: the inclusive jet spectrum, the distribution of the hardest 
(ptji) and second hardest {pt^j'i) jets, (half) the effective mass of the two hardest jets 
{Ht,2), three hardest jets (-ffT.s) and of all jets above 40 GeV (Ht). Also shown on the 
right are the nLO results for the il'-factors. The NLO and nLO (fi) widths correspond 
to the uncertainty due to simultaneous renormalisation and factorisation scale variation 
by a factor of two around a central value fi = Ptji- The f?,LO(-RLs) width shows the 
uncertainty from a variation of i?LS in the range 0.5 < i?LS < 1-5. 



enhancements in the Z+j case. This provides supporting evidence as to their legitimacy. 
A second comment is that the cross sections are large: these observables will be easily 
accessible with a few pb~^ of integrated luminosity at a 7 TeV LHC, allowing for an early 
experimental verification of the large fT-factor for Ht- 

The other observable in the left-hand plot of fig.[IOl Pt,j2, has a very different i^'-factor, 
somewhat below 1. The right-hand plot shows the NLO i^-factors for our full range of 
observables, focusing on a single bin of the left-hand one, from 400 — 500 GeV. The 
pattern that we see here allows us to make some deductions. Firstly, the Ht,2 variable, 
which sums the pts of the two leading jets, is free of large NLO enhancements. It is the 
addition of the third jet in Ht,3 and Ht that brings about the enhancement. A natural 
interpretation is the following: it is common for a third, soft jet to be present due to initial 
state radiation. This third jet shifts the Hj' distribution to slightly larger values, and 
because the distribution falls very steeply, that leads to a non- negligible enhancement. 
This suggests that if, in section [5l we had used effective mass observables with at most 
two objects in the sum, then the nNLO/NLO ratios would have been close to 1. We have 
verified that this is indeed the case. 

The pattern for pt^i and pt^2 in fig- [IH] can also be explained in similar terms: a soft 
ISR emission boosts the hard dijet system, breaking the degeneracy between the p^'s of 
the two hardest jets. It is jet 1 that shifts to larger pt (giving a i^-factor > 1), while jet 



25 





3 




2.5 


II 


2 






o 

_l 


1.5 


o 




g 


1 








0.5 









3 




2.5 


II 


p 






O 
_i 


1.5 


o 




g 


1 








0.5 








1 ' 1 ' 1 ' 1 — 

pp, 7 TeV, anti-kj R=0.7 
NLOJet++, CTEQ6M 




200 400 600 800 
inclusive Pj j [GeV] 



200 400 600 800 
Ptji [GeV] 



200 400 600 800 
Pt,j2 [GeV] 



1 ' 1 ' 1 ' 1 


1 ' 1 ' 1 ' 1 




" (^^rmmihhnnnnnnnnniiim 












rripHHil 1 1 1 1 1 Ti 1 1 1 1 1 1 1 i i 






1,1,1.1 


1,1,1,1 


1,1,1.1 



200 400 600 800 
i Ht,2 [GeV] 



200 400 600 800 
i Ht,3 [GeV] 



3 

2.5 
2 

1.5 
1 

0.5 


3 

2.5 
2 

1.5 
1 

0.5 




200 400 600 800 
i Ht [GeV] 



Figure 11: The nNLO and NLO K-factors relative to the LO predictions, as a function 
of pt (or \Ht, etc.), for the collection of jet observables described in the text. 



2 shifts to lower pt and so it gets a K factor below 1. For the inclusive jet spectrum, and 
for Ht^2, this effect balances out. In addition, final-state radiation from one of the jets 
can cause it to shift to lower pt (becoming the 2nd jet), further reducing the ii'-factor for 
the distribution of Ptj2- 

Of the different variables, it is only the inclusive jet pt and Ht,2 for which there is a 
clear reduction in scale uncertainty in going from LO to NLO. 

Figure [TOT right) also shows the nLO results (including uncertainties both from scale 
variation and from the LoopSim parameter -Rls)- Despite the fact that none of the K- 
factors is parametrically large (except arguably for if^.s and Ht)-, the nLO results are 
remarkably effective at reproducing the pattern of NLO i^-factors, albeit with a small 
systematic shift and generally larger scale uncertainties. One can also verify that, to 
within 10 — 20%, the pt dependence of the NLO ii'-factors is reproduced at nLO. 

Given this success of nLO, and the observed limited convergence of some of the observ- 
ables at NLO, it is interesting to examine what happens at nNLO, where the additional 
3j@NL0 contribution that we require is again obtained using NLOJet++. Results are 
shown in fig. [TTl 

For the inclusive jet spectrum and Ht^2i which already saw large reductions in scale- 
dependence at NLO, the nNLO corrections have essentially no meaningful effect: they 
neither significantly affect the central values, nor reduce the scale uncertainties. For these 
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observables, NLO already converged well, and adding a subset of the NNLO corrections 
without the 2-loop part cannot improve the result. 

For the other effective mass observables, the situation is quite different. With -ffT,3, 
the nNLO result is close to the NLO result and the scale uncertainty is much reduced, 
i.e. this observable seems to come under control at nNLO. In contrast, Ht is subject 
to quite a large further correction, with the central value at nNLO lying outside the 
NLO uncertainty band, and the nNLO uncertainty band (dominated by scale variation) 
only marginally smaller than at NLO. Why is this? Perhaps we are seeing the effect of 
a second ISR emission, which shifts the Ht distribution to even higher values? Given 
that Ht^^ converges and Ht does not, such an explanation is not unattractive. It is also 
consistent with the decrease in i^-factor at low Ht-, where the 40 GeV pt cutoff on the jets 
contributing to the Ht sum will eliminate the ISR enhancement. A definitive conclusion 
would however probably require further study. 

For the remaining two observables, pt,i and pt,2i the nNLO contribution goes in the 
opposite direction from the NLO correction and at low pt it seems that the series fails 
to converge. This is, we believe, closely related to observations of insufficiencies of NLO 
predictions for dijet cross sections in DIS and photoproduction when identical pt cuts are 
imposed on both jets [SI HHl HH HZl HH] (equivalent to integrating the pt2 distribution 
above that cut). The worse convergence at low pt is probably due to the larger fraction of 
subprocesses that involve gluons in the underlying 2 — i- 2 scattering, so that perturbative 
corrections tend to go as (CaOs/tt)" rather than as {CfO-sIt^Y higher pt. 

Considering that we do not have giant NLO i^-factors for the jet processes shown here, 
one may question the validity of the information obtained from the LoopSim procedure. 
An important cross check comes from a comparison with the reference-observable tech- 
nique. Examining fig. [10] (right), one sees two natural reference observables: the inclusive 
jet spectrum and i?T,2! both of which show "perturbative" i^'-factors and small scale 
dependence at NLO. Here we will use (half) the inclusive jet spectrum as the reference 
observable (results with Ht^ would be almost identical). 

Figure fT2] provides a comparison of the LoopSim nNLO results (showing the envelope 
of the scale and i?LS uncertainties) with the reference-observable nNLO results. The com- 
parison is given for all observables except the reference observable itself. The agreement 
between the two methods is striking, with the reference-observable method giving just 
a small shift of the i^-factors relative to the LoopSim results. The shift is identical for 
all the observables, as it has to be: it is simply equal to the difference between the NLO 
and nNLO results for the reference observable. Insofar as we believe the scale depen- 
dence to be representative of the true NLO uncertainty on the inclusive jet spectrum!"] 
the results for the other observables should therefore be good approximations to the full 
NNLO resuhs. 

^-'^In light of the fact that the nNLO uncertainty for the inchisive jet spectrum is larger than the NLO 
uncertainty, it may be that our symmetric scale variation is underestimating somewhat the uncertainties 
present at NLO. To be conservative, it might have been safer to vary the renormalisation and factorisation 
scales independently. 
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Figure 12: Comparison of LoopSim based nNLO results with those from the reference- 
observable method, here using the inclusive jet pt spectrum as the reference observable. 
In the LoopSim results (labelled nNLO), the uncertainty bar spans the envelope of the 
scale and i?LS uncertainties. The results are for the 400 < V/ GeV < 500 bin for each 
observable V, as in fig. [TOl 

7 Conclusions 

Several cases of LHC observables with giant NLO i^-factors have come to light in re- 
cent years. They are characterised by the presence at NLO of new partonic scattering 
topologies that have large enhancements over the LO topologies. In these cases, NLO 
calculations, while important in highlighting the presence of the large i^'-factors, cannot 
on their own provide accurate predictions. 

In this article we have examined how to address this problem by combining NLO re- 
sults for different multiplicities, for example Z-|-j@NLO with Z+2j@NL0. Our main, most 
flexible method, LoopSim, makes use of unitarity to cancel the infrared and coUinear di- 
vergences that appear when one tries, say, to apply Z+2j@NL0 calculations to observables 
that are non-zero starting from Z+l-parton. We referred to the result as Z-f-j@nNLO, 
where the "n" indicates that the highest loop contribution to the NNLO result (the 
two-loop part) has been estimated with LoopSim. 

In introducing a new approximate method for estimating NNLO corrections, sig- 
nificant evidence needs to be provided that the method is meaningful. Firstly, we 
gave reasons why, in cases with giant i^-factors associated with new NLO topologies, 
we expect nNLO results to be a good approximation to NNLO results. As a next 
step, we carried out studies comparing Z/7*@nNL0 (DY) to NNLO predictions for the 
pp ^ Z/'y* + X ^ e~^e~ + X process. In comparing the DY lepton pt nNLO distributions 
to NNLO we found near-perfect agreement in a region of giant i^-factors, pt — |niz > Tz- 
Interestingly, even in the region where the NLO ii'-factor was not large, pt < ^mz, the 
nNLO results provided a significantly better approximation to NNLO than did the plain 
NLO result. This need not always be the case, but is, we believe, connected to the ob- 
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servation that our nLO results reproduced much of the structure seen at NLO (recall, 
Z@nLO means combining Z@LO with Z+j@LO). 

For Z+j production, the first step of our validation procedure was to compare nLO 
and NLO results. All observables with giant /^-factors showed good agreement between 
the two (one with a moderately large /^-factor did not). For those observables, fiNLO 
always appeared to provide extra information: either suggesting a convergence of the 
perturbative series, with reduced scale uncertainties (for ptji), or an indication of sub- 
stantial further higher order corrections (for the effective-mass type observables -f^Tjets 
and if^.tot)- Almost identical results were seen with our alternative "reference-observable" 
estimate of the NNLO contribution. 

The large nNLO corrections that we saw for effective mass observables led us to 
examine a range of effective-mass and jet observables in the simpler context of pure jet 
events (with the expectation that Z+j@NNLO might be similar to 2j@NL0). There we 
saw a significant NLO i^'-factor for all effective mass variables except one, -ffr,2; which 
summed over just the two leading jets. In the Z+j case we had summed over all jets and 
hence it is not surprising that we should have observed substantial nNLO/NLO ratios. 

Even though the observables in the pure jets case did not display giant ii"-factors, the 
pattern of NLO results was remarkably well reproduced at nLO. This encouraged us to 
study nNLO predictions, which provided substantial extra information for several of the 
observables, with the reference-observable method again giving important cross checks. 
Since the cross sections for the jet observables are large, these results could easily be 
tested with early LHC data. 

We close this article with a few lines on the relation between LoopSim and other 
predictive methods. There is a close connection between nLO (or nnLO) and CKKW 
and MLM H5] matching, since they also both provide ways of combining tree-level 
results with different multiplicities. Of course CKKW and MLM matching provide an 
interface with parton showers too, which the LoopSim method does not. On the other 
hand it is significantly easier to include multiple loop orders into the LoopSim method 
than it is within matrix-element /parton-showering matching procedures (though work is 
ongoing in this direction see e.g. [50]). 

An interesting cross-check of the LoopSim method will come with the completion 
of the NNLO calculations for the Z+j and dijet processes. At that point the method 
could also, for example, be used to merge Z@NNLO with Z+j@NNLO, so as to provide 
an nNNLO prediction for quantities like the Drell-Yan lepton pt spectrum. The value 
of the LoopSim method also goes hand-in-hand with progress on 1-loop calculations, 
especially with the prospect of automated NLO calculations now on the horizon (for 
example [SB [S21 [S3]). 

Note that the LoopSim code, which will be made public in due course, can currently 
only deal with hadron-collider processes involving any number of light partons and up 
to one vector boson. It would benefit from further work to appropriately include heavy 
quarks and additional bosons. 
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A Recoil procedure 

In this appendix, we provide further details on how we perform the recoil of an event 
when a particle becomes virtual, including the treatment of the decay products of the Z 
boson. We first examine the simpler case of a particle that makes a loop with the beam, 
then we show how to deal with a particle that makes a loop with another particle. 

A.l A particle recombines with the beam 

Let us assume that particle iq makes a loop with the beam. To balance the event moment, 
we follow the following procedure: 

1. For each particle i ^ iq, store its rapidity Ui. 

2. Perform a separate longitudinal boost on each particle so as to bring its rapidity to 
(i.e. get a purely transverse event). 

3. Compute 

Etot = XI ' (31) 

where Ei is the energy of particle i in the purely transverse event. 

4. Define 

k = {E = Etot,Pt = Pt,io,Pz = 0) , (32) 

and boost all particles into the rest frame of k (so that the total transverse momen- 
tum balances). 

5. Perform a longitudinal boost on each particle so that it recovers its original rapidity 
Vi- 

For the case where two particles, io and ii, are looped with the beam, replace ^ 7^ ^o 
with i 7^ io,ii and in eq. (132|) replace pt^i^ with pt^i^ + Pt,ii, etc. In the case where the Z 
decays, for instance into 2 leptons, the procedure is identical except that we apply to the 
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Figure 13: Case where four gluons are emitted from the same quark. Gluon 1 is the last 
to be clustered with the quark (which roughly corresponds to an early time emission) 
and gluon 4 is the first to be clustered. In the event where gluon 2 makes a loop over the 
quark, we spread the gluon 2's momentum over the quark's momentum and the momenta 
of gluons that were emitted after it, i.e. gluons 3 and 4 (an earlier time emission like gluon 
1 cannot be affected). 



leptons the same longitudinal boosts as for the Z (the rapidity of the leptons is thus not 
necessarily when we apply the transverse boost). This conserves the property that the 
sum of the leptons' momenta is still the Z momentum in the "looped" event. 

The logic of the above procedure is that if we had attempted to apply a transverse 
boost without stages 2 and 5, we would have found that our choice of transverse boost, 
and the corresponding mapping of high-pt particles' momenta, would be affected by the 
presence of energetic particles coUinear to the beam. This would have made the procedure 
coUinear unsafe. 

A. 2 A particle recombines with another particle 

Let us consider the situation depicted in fig. [T3j four gluons are emitted from the same 
quark, but at different angles: 

Oi, > 92, > 9s, > 9,g , (33) 

and gluon 2 becomes virtual. The virtualisation of gluon 2 over the quark cannot have an 
impact on gluon 1, which was emitted earlier in an angular-ordered picture. But it has 
an impact on gluons 3 and 4. More precisely, let the pi be the momenta in the original 
event and the momenta in the event where gluon 2 is virtual. We define 

Pt,tot = Pt,q + Pt,3 + PtA ■ (34) 

and then set the p- as follows: 

p'i=Pi + -^P2 for i = g, 3, 4 , (35a) 

Pt,tot 

p'l = Pi ■ (35b) 

Subsequently each particle's p^ momentum is adjusted such that its mass is (or rriz 
if gluon 3 is a Z boson rather than a gluon), keeping its transverse components p^, Py 
and its rapidity unchanged. This can be easily generalised to any number of particles 
recombining with the same hard one: for each recombined particle i, we spread the looped 
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particle over the hard particle h and over any non-looped emissions from h that are at 
smaller angle (i.e. earlier in the C/A clustering sequence) than i. In eqs. fl34f35ap it is 
always the original particle momenta that are used to determine the pt,i/pt,tot ratio, so 
that the result is independent of the order in which we perform the recombinations. 

This procedure is designed to ensure collinear safety: if, for instance, gluon 4 is 
collinear to the quark in the original event, then it remains collinear in the looped event. 
And if it is gluon 4 (emitted after all the others) that is looped, only the quark momentum 
is rescaled and its direction barely changes, so that angles between the quark and the 
other gluons stays the same. 

In the case where the Z decays into 2 leptons, one applies the following procedure to 
each of the leptons: 

1. Perform a longitudinal boost of the Z boson respectively in the original event and 
the looped event such that it has rapidity in each case. Call the momenta obtained 
pz,o = {Eo,pt^o,0) and pzA = {Ei,pt,i,0) respectively. 

2. Perform a longitudinal boost of the lepton from the original event into the frame 
where the initial Z has rapidity. 

3. Define a purely transverse vector k such that pz,o is transformed to pz,i if it is 
boosted into k's rest frame: 

k=(E, + Eo, - ^*'0)' o) ' (36) 

with 

^" {E, + E,Y ■ ^^^^ 

4. Boost the lepton's momentum into fc's rest frame. 

5. Apply to the lepton the longitudinal boost that brings pz,i to its true rapidity in 
the looped event. 

We are aware of the cumbersome nature of these procedures. A simplification of them 
that retained the relevant coUinear-safety properties would certainly be of interest. 



B The LoopSim method and incoming partons 

Without going into a full proof, we shall here illustrate why the LoopSim method is 
sensible even in the presence of incoming hadrons, by considering what happens at nLO. 
We start with a LO cross section for a process producing n hard objects 

dXadXbd^n '•^^"^^ fi/a{Xa,l4)fj/b{Xb,l4) C{pi, ■ ■ ■ ,Pn) ■ (38) 

For compactness of notation, we have dropped the /ia dependence in the differential tree- 
level partonic cross section daij^n/d^n- We have also not yet specified our choice for the 
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factorisation scale /ip. We assume that daij^n/ d^n contains the necessary constraints to 
relate the incoming partonic momenta to the outgoing momenta. We further integrate 
over the phase-space dxadxbd^n, and include a function C{pi, . . . ,Pn), which is 1 if the 
momenta pass our cuts and otherwise. 

We now imagine that there is some transverse-momentum scale Qq below which no 
radiation is emitted. To 0{as), the PDFs fi/a{xa, fJ''^) can be written in terms of PDFs 
at scale Qq: 

as f^^ dk^ f dz 

ft/aiXa, fJ^l) = fi/a{Xa,Ql) + 7^ / "rf / — Pik{z) fk/a{x / Z , QI) , (39) 

where we sum implicitly over repeated indices. We also define an unregularised splitting 
function Pik{z) such that Pik{z) = Pik{z) — 6{1 — z) f dz'pik{z'), with Pik{z') embodying 
the virtual parts of the splitting function (it is zero for i ^ k). 

Next, we write the LO cross section in terms of a PDF for proton a that has been 
evaluated at scale Qq: 



LO f ^ ^ ddij^niXaPai XbPb) n / 2\ / \ 

- ' dXadXbd^n — fj/b{Xb,fXp)C{pi,...,Pn) 



X 



fi/aiXa,Ql) + ^ I -T^dz ( ^'''^J^ fk/ai^a/z, Qq) - Pik{z) fi/a{Xa, QI)\ 

J Qq t \ / 



(40) 



Note that the first term in round brackets on the second line corresponds to real emission 
of a parton. However that parton is not taken into account in the C{pi, . . . ,Pn) factor. 
Next we examine the structure of the nLO contribution, 

fiT.O TjO . f 1 1 1^ d&ij^n+li^aPay XbPb) e i 1\ r i 1\ 

= + dXadXbd^n+1 — fi/a{Xa, fXp) fj /b{Xb, l^p) 

J Ct^n+l 

X [C{p„...,p,,+^)-C{p'^^...,plf)] , (41) 

where the p\^ . . . p^f represent the momenta when the LoopSim procedure has looped 
Pn+i- In the limit in which Pn+i is coUinear to incoming parton i, with momentum 
Pn+i ~ (1 — z)xaPa, the n+l-partou differential cross section and phase-space simplify 

, ,^ daij^n+l{XaPa,XbPb) _ dakj^nKPa, XbPb) as dz d^f^ 

where x'^ = zxa- By "collinear" we will mean ^ Q where Q is the momentum 

transfer in the hard process. In this limit we also have that pf^ ^ pi (for I <n). So, still 
working within the collinear limit, we can now rewrite eq. fHT]) as 



nLO LO I [jlj dcfkj^niXgPa^XbPb) r / 2\ 
nLU r^cr^^+ I dxJXb d^n fj/b{Xb, /ip) 



a 



^ I -TT^ — Pkt{z)[C{pi,...,Pn+l)-C{pi,...,Pn)]^/a{x'a/z,fJ.l). (43) 
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Next, we exchange i k, replace x'^^ — )■ Xa and change the scale in fijaix'^j z, /ip) to be 
Qq, which is allowed because it corresponds to an 0{o?^ change (while here we consider 
only 0{as)): 



Cyf'^^O^^^ I dxjx,d<l>r. da,,^4x^p^^,X,p,) ^^.^^(^^^^2) 



^ -JT Pik{z)[C{pi,. . . ,Pn+l) - C{pi,. . .,Pr,)]fi/a{Xa/z,Qo) . (44) 

^TT Jq2 ftj^^^i Z 

Note now that if we take /ip ~ in eq. fHOj) . then the second term in square brackets in 
eq. (jH]) cancels the first term in round brackets in the second line of eq. fHOj) . In other 
words for initial-state radiation, the action of LoopSim is not so much to provide virtual 
corrections as to cancel the real-emission terms already included implicitly through the 
PDFs in the leading order cross section. In contrast, the true virtual terms are already 
included through the PDFs themselves, i.e. through the second term in round brackets 
in eq. f l40|l . 

As an example, consider pp — Z. At nLO we will have events such as gq — )■ Zq, where 
the outgoing quark comes from coUinear initial-state splitting g — )• qq, with an underlying 
hard subprocess qq — )■ Z. From these events LoopSim will generate a configuration in 
which the outgoing quark is "looped" . This will come in with a PDF weight that is the 
product of a gluon distribution and a quark distribution, so it appears that we have a 
(negative) gq ^ Z contribution, which would be unphysical. However in the LO cross 
section with a factorisation scale fi^ ~ Q, when we write qq — )■ Z, part of the q PDF 
comes from g ^ qq splitting. If we were just to add the real gq — ?■ Zq diagram to 
the LO cross section alone, then in the collinear limit we would be double counting the 
part already included in the PDF. With the negative "(^g — )■ Z" LoopSim contribution, 
what happens is that we simply remove the q PDF component, generated from g ^ qq 
splitting, that was implicitly included at LO with an incorrect final state (i.e. lacking an 
outgoing quark), since we are now putting it in with the correct final state through the 
real gq — > Zq diagram. 

Note that we have not yet worked out the full extension of this discussion to higher 
orders. The details would depend on the precise higher orders that we have in mind, for 
example nriLO versus fiNLO. However, regardless of these details, the fundamentally 
unitary nature of the LoopSim procedure is important in ensuring that the simulated 
"loops" simply bring about an overall consistent set of final states while maintaining the 
total cross section as calculated with a sensible factorisation scale choice. 



C Secondary emitters in LoopSim 

In section 12.1.21 we discussed the special treatment needed for "secondary emitters", i.e. 
non-Born particles that have emitted something. In our procedure, secondary emitters 
do not get looped: when particle j makes a loop over i, this is justified by the collinear 
enhancement of the matrix element due to i and j being close in angle. But the emission 
of the same j from the configuration where i is virtual does not have such a collinear 
enhancement, so one must not take it into account. Another way to understand it is to 
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consider the example of 2-gluon emission from a qq dipole. The squared matrix element 
for the emission of 2 real energy-ordered {Ei ^ E2) gluons, gi, g2, can be expressed as 

M{h, k2) = {AnasfiC^pWi + C^C^W^s) , (45) 

with 

^ {Pq.ki){ki.pq) \{Pq.k2){k2-kl) {Pq-k2){k2-kl) {Pq-k2){k2-Pq) J 

Since the Wi term diverges when g2 is collinear to g or g (unlike the W2 term), it becomes 
relevant when g2 is considered to have been emitted from g or g independently of gi . The 
W2 term diverges when g2 is collinear to gi (unlike Wi), so it becomes relevant when g2 
is considered to have been emitted from gi. This is depicted in fig. [TH which also shows 
the virtual corrections (cf. [55] )• One notices that the W2 term only appears when gi is 
real. The diagrams where gi is virtual are taken into account when g2 is emitted from g 
or g. Therefore, gi cannot become virtual when (72 makes a loop over it. 
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Figure 14: Schematic depiction of the matrix elements for two gluons emitted (or virtual) 
from a gg dipole: (a) gluons 1 and 2 real; (b) gluon 1 real and 2 virtual; (c) gluon 1 
virtual and 2 real; (d) gluons 1 and 2 virtual. In each case, when needed, we use the 
decomposition into Wi and W2 pieces to separate what can be seen as the emission of 
gluon 2 from gluon 1, and what can be seen as the emission of gluon 2 directly from gg. 
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